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We  show  that  the  low-lying  spin  states  of  two  electrons  in  a  semiconductor  quantum  dot  can  be  strongly 
mixed  by  electron-electron  asymmetric  exchange.  This  mixing  is  generated  by  the  coupling  of  the  electron  spin 
to  its  orbital  motion  and  to  the  relative  orbital  motion  of  the  two  electrons.  The  asymmetric  exchange  can  be 
as  large  as  50%  of  the  isotropic  exchange,  even  for  cylindrical  dots.  The  resulting  mixing  can  contribute  to 
understanding  spin  dynamics  in  dots,  such  as  recent  observations  of  light  polarization  reversal. 
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An  electron  spin  in  a  semiconductor  quantum  dot  (QD)  is 
an  attractive  qubit  for  quantum  computing:1  the  spin  in  the 
ground  orbital  state  can  have  a  long  coherence  time;2  a  single 
qubit  can  be  initialized  or  read  optically  by  transient 
electron-hole  pair  excitation  giving  a  negative  trion  A-;3,4 
and  the  manipulation  of  the  spin  exchange  between  spins  can 
be  the  basis  for  two-qubit  gates.1  The  correlations  between 
spins  control  coherence  in  gates.  The  dominant  interaction 
between  two  electrons  ( e-e )  is  the  symmetric  exchange 
Jsy§2,  which  conserves  the  total  spin  S.  Additional  spin- 
asymmetric  e-e  interactions  (asymmetric  exchange)  do  not 
conserve  S  and  can  cause  decoherence. 

Examples  of  recent  experiments  on  spin  dynamics  are 
those  involving  optical  polarization  reversal.3,4  They  involve 
spin  flipping  due  to  electron-hole  ( e-li )  exchange  in  QDs 
with  lateral  asymmetry.  It  has  been  noted  that  those  experi¬ 
ments  require  strong  spin  mixing,  inconsistent  with  e-li  ex¬ 
change  alone.4 

Spin-orbit  (SO)  interactions  play  a  key  role  in  the  mixing 
of  spin  states.  They  arise  from  effective  magnetic  fields  of 
the  orbital  motion  of  electrons.5  Electrons  in  QD  ground 
states  with  dominant  s  components  have  small  orbital  angu¬ 
lar  momentum  and  small  SO  coupling.  A  number  of  experi¬ 
ments  involve  electrons  in  excited  states.  Linear  combina¬ 
tions  of  nearly  degenerate  excited  states  in  a  plane  (e.g.,  px- 
and  /7,,-like)  can  give  rise  to  two-dimensional  (2D)  orbital 
motion  with  an  effective  magnetic  field  perpendicular  to  the 
plane,  and  thus  to  large  SO  coupling.  This  is  analogous  to  the 

L  S  coupling  in  atoms.  Thus,  symmetric  QDs  (e.g.,  cylindri¬ 
cal)  can  have  significant  SO  effects. 

There  are  three  sources  of  SO  coupling  that  lead  to  mix¬ 
ing  of  spin  states.  The  largest  two  contributions  arise  from 
the  kp  mixing  of  the  conduction  and  valence  bands  and  can 
be  described  in  the  effective  mass  approximation.6  We  derive 
them  by  treating  the  potentials  from  the  structure  and  the  e-e 
Coulomb  repulsion  on  the  same  footing  with  k  p  terms  using 
the  Kane  model.7  We  have  in  mind  QDs  with  a  strong  con¬ 
finement  in  a  single  state  £(z)  in  the  vertical  direction  e:  and 
weaker  confinement  in  the  transverse  directions,  which  gives 
the  single-particle  electron  states  </>,■(/")  =  £(z)  <Pi(p)  • 

A  single-electron  contribution  to  the  SO  coupling,  hv, 
arises  from  the  2D  structure  potential  V(r):8 

hV-s  =  JVs[d7V(pL  X  s4)  +  ( dpV  X pL)sz]ez.  (1) 

pz  is  not  present  due  to  the  strong  vertical  confinement  [for  a 
single  state  £(z),  (£|/r|  £)  =  0].  The  first  term  is  the  usual 


Rashba  coupling  yv(ez  X/? 1 )  from  the  asymmetry  in  the 
growth  direction,9  where  yv =  yvs{£\  dzV\  f).10  The  second 
term  is  important  between  (almost)  degenerate  excited  states, 
e.g.,  px-  and  p,,-like,  where  it  gives  the  dominant  SO  cou¬ 
pling  independent  of  structure  or  bulk  inversion  asymmetry; 
it  is  negligible  in  the  ground  state,  whose  main  component  is 
inversion  symmetric  (s-like). 

The  second  contribution,  hc,  is  from  the  interaction  of 
each  spin  with  the  other  electron’s  orbital  motion.  We  obtain 
it  using  a  two-particle  k  p  approach  for  electrons  interacting 
through  the  Coulomb  potential 11  Uc(rr)  =  e2/ kt,.  (k  is  the  di¬ 
electric  constant;  rr=r\-r2 ): 

hk  ■  4  =  (- 1 )kys(yruc  X  pk)  ■  sk,  (2) 

^  T  T 

where  k=  1,2.  h  is  analogous  to  the  Pauli  SO  interaction, 
while  hc  is  analogous  to  the  Breit-Pauli  spin-relative-orbit 
coupling.5  These  couplings  in  vacuum  or  in  atoms  are  rela- 
tivistically  small  due  to  the  large  energy  gap  2 m0c2  between 
electron  and  positron  bands,  whereas  the  present  gap  Eg  is 
smaller  giving  larger  SO  couplings. 

A  smaller  contribution,  hB,  comes  from  the  Dresselhaus 
coupling  due  to  the  lack  of  bulk  inversion  symmetry.12  It 
arises  from  the  mixing  of  the  conduction  band  with  the 
remote  upper  bands.  In  QDs  with  strong  vertical  confinement 
[(/r2)»(/?42)],  it  reduces  to  hB-L—  yB{px, -/?,,),  where  y® 

=  Mp2Mn 

We  use  a  model13  of  QDs  like  those  from  self-assembled 
growth  along  the  crystal  axis  [001].  The  lateral  potential 
V(p)  contains  Vs  symmetric  for  the  inversion  p— >-p,  and  it 
may  also  contain  an  inversion-asymmetric  part  Va.  The  lat¬ 
eral  parameters  Dxy  of  the  potential  provide  a  measure  for 
the  QD  size.  A  nonzero  asymmetry  Va=Vax+Vay  (Vflt  odd  in 
x,  and  Vay  odd  in  y)  implies  a  nonzero  average  lateral  electric 
field  £x  y.  The  inversion  asymmetry  is  parametrized  in  Vax  by 
Ex,  and  in  Vay  by  Ey.  We  can  write  £xy°^Exy/Dxy.13 

First,  we  consider  the  orbital  eigenstates  of  the  two- 
particle  Hamiltonian  H0  containing  the  Coulomb  interaction 
U c  but  not  the  SO  couplings.  //0  has  an  inversion-symmetric 
part  (pf2+p22)/2m  +  Vs(Pi)  +  Vs(Pi)  +  Uc{r,)  +  yc  S(r ,.) ,  and 
an  inversion-asymmetric  part  Va(p\)  +  Vfl(p2).  From  their 
permutation  symmetry,  the  eigenstates  are  separated  into 
triplets  {I1,}  (asymmetric)  and  singlets  (5)}  (symmetric).  We 
build  a  basis  for  each  group  from  products  of  single-particle 
wave  functions.  We  use  a  large  set  of  harmonic  oscillator 
wave  functions  obtained  from  the  average  curvatures  of  the 
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FIG.  1.  (Color  online)  QDs  with  diameters  Dx=Dy=D.  (a)  Two- 
electron  energy  levels  in  QDs  with  D=  19  nm  with  one  plane  of 
symmetry  along  ex  (Ey= 0)  vs  the  lateral  asymmetry  parameter  Ex 
(the  average  field  is  £X=\0.\EX  meV/nm).  (b)  The  asymmetric  ex¬ 
change  lycfjg  |  =fi~l  and  its  components  \p.z{3\  (from  the  SO 

coupling  [Eq.  (1)]),  \p-\3\  (from  the  spin-relative-orbit  coupling 
[Eq.  (2)])  vs  the  QD  size  D.  (c)  |/r/j3|  vs  Ex  (Ey= 0)  for  several  QDs. 


(c) 


potential  and  grouped  in  four  subspaces  {s},  {x},  {y},  and  {d}, 
by  their  orbital  symmetry:  s  symmetry  (even  in  x  and  v),  x 
(odd  in  x),  y  (odd  in  y),  d  (odd  in  x  and  y).  The  triplet  basis 
{t„,}  contains  asymmetric  combinations  such  as  tsx=(sx 
-xs)/\2,  and  has  four  independent  subspaces,  of  symmetry 
s,  x,  y,  and  d:  {rss, ,  txx,,  ryy, ,  rdd,},  {rsx,ryd},  {rsy,Txd}, 
{rxy,  Tsd}.  The  orbital  singlet  basis  {cr„}  is  ^iven  by  symmetric 
combinations  such  as  crss/  =  (ss'  +  s's)/ V2  (sAs')  and  trss 
=  ss,  and  it  has  four  subspaces  {crss,,mxx/,cryy,,o-dd,}, 
{crsx,cryd},  {<xsy,crxd},  and  {crxy, <Tsd}.  These  subspaces  are  not 
mixed  by  the  inversion-symmetric  part  of  H0.  The  inversion- 
asymmetric  part  of  770  couples  the  singlet  subspaces  among 
themselves  by  terms  linear  in  Exy  and  the  triplet  subspaces 
among  themselves.  The  eigenstates  of  H0  have  several  com¬ 
ponents  from  different  subspaces: 

T,=T]  +  ExT\  +  ExEyTy  +  Eyr\, 

S3  =  Sl  +  EyS’l  +  ExEyS*  +  ExSi  (3) 

and  similarly  for  72  and  S2-  Here  7j  (72)  labels  the  lowest 
triplet  with  a  dominant  x  (y)  component  and  S2  (S3)  labels 
the  lowest  singlet  with  a  dominant  x  (y)  component.  T\  is  the 
projection  of  7j  on  the  x-symmetry  triplet  subspace,  etc.  The 
lowest  states  are  shown  in  Fig.  1(a)  where  Ey=0.  Higher 
states  not  shown  are  T3  (74),  which  are  the  lowest  d-  (,v-) 
dominant  triplets,  and  S4  (the  lowest  <7- dominant  singlet). 
The  isotropic  part  of  the  exchange  for  triplet  7)  (energy  ej) 
and  singlet  Sj  (energy  ej)  is  given  by  7,y=2(e[-  ej)/7i2.  In  this 
work  we  choose  the  energy  splitting  between  the  electron 
ground  and  excited  states  to  be  in  the  range  20-45  meV;  this 
gives  an  exchange  splitting  (713  between  7)  and  .S’3 )  of  the 
order  5-10  meV,  in  the  range  of  experiments.4 

Next,  the  triplet-singlet  mixing  comes  from  the  SO  terms 


hv,  hc,  hB  added  to  770.  These  give  a  Hamiltonian  composed 
of  a  spin-symmetric  part  77,  that  conserves  the  total  spin  S 
=  V|  +s2  and  a  spin-antisymmetric  part  Ha : 

Hs  =  770  +  ^(h ,  +  h2  +  ysdpUc  X p))  ■  S, 

Ha  =  i  -  h2  +  2Js')„Uc  X  p))  ■  (s  i  -  s2),  (4) 

where  hk=hl+hB,  pr=p\-p2,  and  pc=(p\+p2)/ 2.  Ha  can  be 
written  as 

Ha  =  2  Pij  ■  (5,  -  s2) I T,)(Sj\  +  H  .  c  .  ,  (5) 

ij 

where  pii=(Ti\hl  + ysdp  UcXp^\Sj)  gives  the  asymmetric 

exchange.  States  of  different  total  spin  S  are  coupled  via  the 
operator  Sl-s2  equivalent  to  the  Dzyaloshinskii-Morya  form 
2H h(s1  Xs2).14  We  can  write 

P  •  (f,  -  s2)  =  PV\  ~  ®  +  PL  ■  ~  it)  ■  (6) 

Pz  conserves  the  total  spin  projection  Sz,  i.e.,  it  mixes  sin¬ 
glets  with  triplets  with  S:= 0  (“longitudinal  mixing”).  This  is 
equivalent  to  a  precession  of  the  total  spin  around  ez  (A Sz 
=  0).  pL  mixes  states  with  different  spin  projections  (|A5':| 
A  0)  (“transverse  mixing”).  The  degree  of  triplet-singlet  mix¬ 
ing  is  given  by  the  ratio  of  the  asymmetric  to  the  symmetric 
exchange:  p.)=/i  'PJJfj.  PtjJ J ij- 

We  group  the  operators  giving  /7,y  in  Eq.  (5)  into  an  axial 

vector  operator  A  =  Ae_  and  two  polar  vector  operators  P 
=  Pez,  R=Rl. 

A  =2yvs(dpVs  X  dp)  -  2  ys{dpUc  X  dp), 

P  =  2yvs(dpyaXdp), 

R  =  -2 yv{ez  X  dp)  +  2 yB(exdXi  -  eydy) .  (7) 

A  and  P  include  the  vertical  magnetic  field  from  the  2D 
motion  in  the  nearly  degenerate  excited  states,  and  they  gen¬ 
erate  pz.  R  arises  from  the  Rashba  and  Dresselhaus  terms, 
and  it  generates  pL. 

Table  I  gives  the  matrix  elements  between  T,  (7=1, 4)  and 
Sj  (/'=  1 ,4)  [Eq.  (3)]  from  the  spin  mixing  operator  A  in  Eq. 
(7).  The  states  are  characterized  by  the  symmetry  of  their 
dominant  wave  function  components,  e.g.,  S2  ~x  symmetry. 
Table  II  gives  corresponding  results  from  P.  The  terms  in 
small  boxes  on  the  second  diagonal  in  Table  I  are  dominant 
and  are  independent  of  lateral  asymmetries.  All  the  other 
terms  in  Tables  I  and  II  are  nonzero  only  for  cases  of  lateral 
asymmetry.  The  central  2X2  block  highlighted  is  of  interest 
for  the  dynamics  of  X~  in  the  “p”  shell.  ’1 

The  matrix  elements  in  Tables  I  and  II  can  be  understood 
by  writing  the  operators  in  the  basis  {rm, cr„}:  A+P 
=1,mtn(Amn+Pmn)ez\ rm)((rn\  +H.c.  The  matrix  elements  A")3 
and  P'j)''  in  the  tables  are  given  by  sums  of  matrix  elements 
Amn  and  Pmn,  respectively,  with  the  same  symmetry.  From 
Eq.  (7),  it  is  seen  that  Amn  is  nonzero  only  for  |Tm)(cr„|  odd 

both  in  x  and  in  y.  Thus  A  can  produce  longitudinal  mixing 
P]j  between  two-electron  eigenstates  7)  and  .Sj  if  one  of  these 
contains  an  x-  (,v-)  symmetry  component,  and  the  other  has  a 
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TABLE  I.  The  part  of  the  longitudinal  coupling  fijj  from  A  [Eq.  (7)].  The  matrix  elements  A f^=(Tf\A\ sf)  are  between  components  of 
definite  symmetries  Tf,  S ?  of  the  orbital  wave  functions  Th  Sj  [Eq.  (3)]. 


S 1  (a;  3-symmetry) 

S2  f~  x-symmetry ) 

S-j  («  y-symmetry) 

S 4  (a:  d-symmetry) 

Ti 

(w  s) 

ExEv[A^+A^l+AiU^] 

Ex  [ A*J  +  A% + EI<.A**+a£ ) 

ah  + 

El  +  E2y(  AH  p  El  AH) 

r, 

(«  X ) 

EzEu\A*Z+A'(Z+Ati+A,ti] 

Ki  + 

EiA^+EKA^  +  ElA^) 

EAAllPAli+EllAH+AlDi 

t2 

(»  y ) 

Bm  [A*f  +Mi + El  (A?l  +  A% )) 

Kz  + 

ElAl‘2+El<,A£+ElA%) 

ExEv[All+A^A^+Ai  I) 

Ey\AX+AlUEl(A'‘l^A™)] 

(«  d) 

All  + 

E*A»*+El(Ax3*+ElA*31) 

EAA&+A'£+El(A£+Ali)} 

Ey  t  AH  +A%+ El  (Ag + Ag)) 

ExEv[All+A»l+A£pAll\ 

y-  id-)  symmetry  part  (Table  I).  Pmn  is  nonzero  only  for 
nonzero  QD  asymmetries  (Va  A  0)  and  for  |rmXtrn|  odd  either 

only  in  x  or  only  in  y.  Thus,  P  contributes  to  the  longitudinal 
spin  mixing  /3(;  between  T,  and  Sj  if  one  of  them  has  an  s  or 
d  component  and  the  other  has  an  x  or  v  component  (Table 
II).  R  can  be  written  as  R=1,m„R^nn\  rm)(crn|  +H.c.  Results 
for  the  matrix  elements  of  R^m  are  not  given  explicitly  here. 
They  require  QD  lateral  asymmetry  and  are  nonzero  for 
I  Tm)(an\  °dd  in  one  of  x  or  y.  They  can  give  transverse  spin 
mixing  of  states  with  different  z  spin  projections. 

We  now  consider  QDs  with  several  symmetries  and  the 
longitudinal  spin  mixing  fiz  from  them.  This  longitudinal 
mixing  does  not  have  contributions  from  the  Dresselhaus  and 
Rashba  couplings. 

(i)  QDs  with  lateral  inversion  symmetry  (EX=E  =  0).  Ex¬ 
amples  are  shown  in  Fig.  1(b)  and  by  the  Ex=  0  points  in 
Figs.  1(a)  and  1(c)  and  Fig.  2(a).  In  such  QDs,  the  two- 
electron  states  ThSj  [Eq.  (3)]  have  well-defined  symmetries. 
The  spin  mixing  is  due  only  to  A  on  the  second  diagonal  (in 
small  boxes)  in  Table  I.  “Pure”  states  of  x  (y)  symmetry  such 
as  Tl  (7  2)  couple  only  to  “pure”  states  of  y  (r)  symmetry 


such  as  S3  (S2).  The  first-order  longitudinal  spin  mixing  of  7) 
(7’2)  is  given  by  the  coupling  S3  (S2),  because  this  is  the 
closest  in  energy.  7’3  (the  lowest  cZ-symmetry  triplet)  couples 
by  A  to  s-symmetry  singlets  like  5j.  T4  (the  lowest 
i-symmetry  triplet)  couples  by  A  to  rZ-symmetry  singlets 
such  as  S4. 

From  Figs.  1(b)  and  2(a)  (at  Ex=  0),  the  asymmetric  ex¬ 
change  can  be  a  substantial  fraction  of  the  symmetric  ex¬ 
change  (up  to  —50%).  From  Fig.  1(b),  the  asymmetric  ex¬ 
change  is  smaller  for  larger  QDs,  which  results  from  larger 
orbits  giving  smaller  effective  magnetic  fields  in  the  SO  cou¬ 
pling.  In  this  case  P22-  -/3j3  because  of  degeneracy.  The  or¬ 
bital  momentum  L,  eigenstates  (S2±iS3)/  f2  are  strongly 
coupled  to  (7)  ±iT2)l  \  2  and  L.  is  conserved.  From  Fig.  2(a), 
the  asymmetric  exchange  decreases  as  the  degeneracy  of  the 
first  two  excited  states  is  removed  when  Dx  A  I), .  In  this  case 
L,  is  not  conserved.  The  stronger  confinement  along  ey  (Dx 
>7?v)  leads  to  J\2>J22,  and  thus  to  |/rj3|  <|/r;2|. 

(ii)  QDs  with  a  single  vertical  plane  of  reflection  (Exdz0 
and  Ey= 0).  This  gives  more  nonzero  matrix  elements  in 
Tables  I  and  II,  e.g.,  now  74  is  mixed  with  .S'3  as  well  as  with 


TABLE  II.  The  part  of  the  longitudinal  coupling  fromP  [Eq.  (7)].  P“P=(Tf\P\S^)  are  between  wave  function  components  of  definite 
symmetry  [Eq.  (3)]. 
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FIG.  2.  (Color  online)  (a)  ptz  in  QDs  with  Dx= 23  nm  and  Dy 
=  15  nm  compared  to  the  longitudinal  coupling  in  QDs  with  Dx 
=DX  =  19  nm.  (b)  Mixing  of  states  with  different  spin  projections  Sz 
in  QDs  with  Dx=Dy  and  with  a  plane  of  symmetry  (Ev  =  0,  Ex  A  0): 
p\2  for  the  mixing  of  7j  and  S2- 

S4.  This  case  corresponds  to  Figs.  1(a),  1(c),  and  2(a).  /i]3  for 
the  lowest  triplet  decreases  with  increasing  Ex.  For  these 
cases,  the  terms  proportional  to  Ex  and  Ex  in  Table  I  and  also 
the  terms  from  P=2yv(dpVax  X  dP]Vmnz>:Ex  from  Table  II  are 
nonzero,  and  they  tend  to  cancel  partially  the  larger  terms  in 
the  boxes  in  Table  I.  For  some  triplet-singlet  pairs,  such  as  S3 
and  7],  the  symmetric  exchange  becomes  larger  and  thus 
their  mixing  decreases.  Other  singlet-triplet  pairs  can  be  de¬ 
generate,  such  as  T2  and  S2  in  Fig.  1(a)  at  Ex~  1.5;  then  a 
nonzero  f£,2  leads  to  strong  singlet-triplet  mixing  (Fig.  2(a)). 
For  this  case,  Lz  is  not  conserved.  Triplets  with  (L:)  ~  ±fi  can 
be  coupled  to  singlets  that  have  (L,)~  +h. 

(Hi)  QDs  with  no  vertical  plane  of  reflection  (Exf=0  and 
Ex  A  0).  Then  all  states  in  Tables  I  and  II  are  mixed,  and  the 
longitudinal  spin  mixing  can  be  larger  than  in  previous  cases. 

In  addition  to  the  longitudinal  spin  mixing  above,  there  is 
also  mixing  that  changes  the  spin  projection  Sz  (transverse 
mixing  /jl  1 ).  It  is  exclusively  from  the  Dresselhaus  and 
Rashba  couplings,  which  give  R  in  Eq.  (7).  For  lateral  inver¬ 
sion  symmetry,  R  mixes  states  which  typically  differ  by  the 
single -particle  energy  splitting,  e.g.,  7j  with  5j  and  .S'4,  etc. 
Then  the  mixing  from  R  is  small,  due  to  large  Jn  and  /14. 
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For  QDs  with  only  one  vertical  plane  of  reflection,  R  mixes 
7’ |  with  .S' i  or  S3,  which  are  closer  in  energy  and  therefore 
give  larger  mixing.  We  show  in  Fig.  2(b)  this  transverse  spin 
mixing  for  T 3  and  S2.  It  occurs  only  for  nonzero  asymmetric 
potential  (Exf=0).  It  is  generally  smaller  than  the  longitudi¬ 
nal  spin  mixing  but  can  become  appreciable  for  large  asym¬ 
metries,  and  it  is  larger  in  smaller  QDs. 

We  have  shown  that  when  the  splitting  between  the  elec¬ 
tron  p-states  in  a  QD  is  small  there  can  be  strong  mixing 
between  electron  excited  singlets  and  triplets.  This  mixing 
can  be  important  in  optical  manipulations  of  spins  in  QDs 
and  can  lead  to  dephasing  and  loss  of  fidelity  in  gates.  For 
example,  these  results  can  help  in  interpretations  of  light 
polarization  reversal  experiments.  These  interpretations  in¬ 
volve  electron  singlet-triplet  mixing.  A  contribution  to  the 
latter  comes  from  e-h  axially-asymmetric  exchange  in  later¬ 
ally  asymmetric  QDs.  The  e-e  asymmetric  exchange  pre¬ 
sented  here  provides  an  additional  contribution,  which  is  im¬ 
portant  in  particular  for  laterally  symmetric  QDs.  The 
present  results  also  suggest  an  additional  process  in  such 
experiments.  Typically  in  optical  pumping,  the  orbital  angu¬ 
lar  momentum  from  light  is  stored  in  the  hole  motion  and  the 
electron  is  in  an  s  state.  For  states  with  an  excited  electron, 
the  electrons  can  carry  the  orbital  angular  momentum.  The 
e-e  asymmetric  exchange  can  mix  two-electron  states  that 
differ  in  their  orbital  angular  momentum  leading  to  emission 
of  light  with  reversed  polarization. 

Finally,  in  gated  QDs,  which  are  typically  larger,  the  mix¬ 
ing  discussed  here  between  triplets  and  excited  singlets  is 
smaller  [Fig.  1(b)].  Nevertheless,  in  situations  like  spin  trans¬ 
port  with  a  bias  through  gated  QDs,15  the  lowest  singlet  can 
be  brought  close  to  the  ground  state  singlet  and  the  mixing 
between  them  can  become  important. 

We  are  grateful  for  discussions  with  Y.  Lyanda-Geller,  D. 
Gammon,  B.  V.  Shanabrook,  and  M.  E.  Ware.  This  work  was 
supported  in  part  by  the  ONR  and  by  DARPA. 
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